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ABSTRACT 

We investigate the formation of molecular hydrogen (H2) in a primordial 
H+He gas cloud irradiated by a power-law UV flux. We find that at high 
densities ( ^ f cm~"^) and low temperatures ( ^ fO'* K), the background radiation 
enhances the formation of H2 and results in molecular cooling dominating 
over photoionization heating. This process could accelerate the collapse and 
fragmentation of dense objects at high redshifts. 

Subject headings: cosmology:theory-early universe-molecular 
processes-galaxies:formation 

submitted to ApJ 
1. Introduction 

It has long been hypothesized, based on studies of QSO absorption spectra (Gunn 
& Peterson 1965), that the intergalactic gas is photoionized by an intense UV background. 
Indeed, a background fiux level ~ I0~^^^^ erg cm~^ Hz"^ sr~^ at the Lyman limit has been 
inferred from the proximity effect in the vicinity QSOs at redshifts 2; 2 — 4 (see Bechtold 
1994, and references therein). This UV background could have resulted from either QSO 
or stellar emission (e.g., Miralda-Escude & Ostriker 1990; Meiksin & Madau 1993; Shapiro, 
Giroux, & Babul 1994; Giroux, Fardal, & ShuU 1995). Aside from its impact on the ionization 
state of the universe, the UV background heats the intergalactic gas and therefore raises its 
Jeans mass. Couchman & Rees (1986) realized the important negative feedback that this 
process might have on the formation of low-mass objects at high redshift. More recently, 
Efstathiou (1992) suggested that photoionization heating of the intergalactic medium to 
temperatures ^ 10"* K could have suppressed the formation of dwarf galaxies (~ 10^ Mq) 
at moderate redshifts, and Babul & Rees (1992) have argued that the formation epoch of 
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faint blue galaxies was indeed delayed by this process until z 1 (see also Babul & Ferguson 
1995). The negative feedback that photoionization heating has on structure formation on 
small mass scales was verified recently by a number of numerical simulations (Hernquist et 
al. 1995; Quinn, Katz, & Efstathiou 1995; Steinmetz 1995; Thoul & Weinberg 1995). 

Naively, one would expect that a relatively cold gas cloud ( ^ 10"* K) could only heat 
due to its interaction with the ionizing background. However, this expectation ignores the 
formation of molecular hydrogen (H2). In this Letter we show that the increased abundance of 
free electrons due to photoionization tends to catalyze the formation of molecular hydrogen, 
which at sufficiently high densities results in a net cooling of the gas. To demonstrate this 
effect we consider a stationary homogeneous gas cloud which is being illuminated by an 
ionizing fiux. §2 describes our method of calculation and §3 presents our results. Finally, we 
briefiy sketch the potential cosmological implications from this effect in §4. 



2. Method of Calculation 



We consider a homogeneous gas cloud of primordial composition at a density p=ijm^n 
and temperature T, made up of the nine species H, H~, H"*", He, He"*", He"'""'", H2, H2', and 
e~. The chemical composition of the gas changes in time according to the rate equations 

dn- 9 9 9 

-17 = aimihmnin^ + ^ kik^j, (1) 

1=1 m=l j = l 

where rii is the number density of the z*^ species, ki^^ denote reaction rate coefficients, k^^ 
denote photo-ionization or photo-dissociation integrals, and the factors aimi, 6ni = 0,ibl, or 
±2, depending on the reaction. The full list of reactions and rate coefficients that we adopt 
are given in Appendix B of Haiman, Thoul, & Loeb 1995. 

We assume that the gas is illuminated by a uniform background radiation with a 
power-law flux extending up to 40keV, 

= J21 (—) X 10-2^ erg s"^ cm-2 Hz"^ sr"^ hp < 40keV. (2) 

Unless stated otherwise, in our calculations we use the values J21 = 0.1 and a = 0.7. This 
normalization of J21 is preferred by recent simulations of the Lya forest (Hernquist et al. 
1995; Miralda-Fscude et al. 1995) and is also compatible with the expected emission from 
the known QSO population (Meiksin & Madau 1993). The power-law index a = 0.7 is 
in accord with quasar spectra, and is also consistent with the flt to the X-ray background 
spectrum around up to 40 keV (Fabian & Barcons 1992). 
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The incident spectrum in equation (2) characterizes the background flux before it 
enters the gas cloud. However, if the gas cloud is optically thick at the Lyman-limit, the 
shape of the spectrum will be modifled due to H-photoionization as the radiation penetrates 
into the cloud. To mimic the effect of radiative transfer, we modify the spectrum to the 
form, 

= (—) exp{-Nn(Jn(iy)} x fO'^^ erg s"^ cm'^ Hz"^ st-\ (3) 

where A^h is the column density (in units cm~^) of neutral hydrogen through which the 
incident flux is processed. For simplicity, we assume that the incident radiation is effectively 
processed through a uniform slab of gas of length £ and deflne 

Here /h is the neutral hydrogen fraction by number, fi is the mean molecular weight, n 
is the total proton number density, and i'jeans is the Jeans length for the assumed density, 
temperature and composition of the slab. 

In addressing the question whether the gas cools or heats due to the radiation 
background at different densities or temperatures, we assume chemical equilibrium. This 
assumption makes our conclusions independent of the details of hydrodynamic and chemical 
history of the gas. (We examine below, as a separate question, the time-scale required 
for establishing the assumed chemical equilibrium). The chemical equilibrium abundances 
of the various species are solved by integrating equations (I) until the right-hand-side 
vanishes. Using these abundances, we then calculate the total cooling and heating functions, 
Acooi and Aheat! these quantities are both positive and are expressed in units erg s~^ cm~"^. 
The interaction between the radiation background and the gas includes the following 
processes: photoionization and dissociation, Compton scattering, coUisional ionization 
and excitation (including H2), recombination, and Bremsstrahlung emission. Numerical 
expressions for these processes are listed in Appendix A of Haiman et al. 1995. In addition, 
for the present calculations we included Compton scattering on the power-law spectrum 
(cf. eq. 24.42 in Peebles 1993). Compton cooling by the cosmic background radiation is 
negligible at 2; ^ 5. 

3. Results and Discussion 

Under chemical equilibrium and the flxed flux parameters (J21 = O.I, a = 0.7), the 
three parameters that deflne the state of the gas are its density p and temperature T, and 
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the column density of the shielding slab, A^h, in equation (3). At each density, temperature, 
and column density, the slab width £ can be expressed in units of the Jeans length using 
equation (4). We first focus on the sign of Atot = Aheat — Acooi as a function of these 
parameters. For reference, the equilibrium H2-fractions as a function of density, temperature, 
and column density A^h are shown in Table 1. Figure 1 shows contours of Atot = in the 
density-temperature plane for various values of neutral-H column densities (A^h)- At high 
temperatures (T ^ 10"* K), H2 has no effect, and cooling always dominates through H- 
excitation. The prominent feature introduced into this figure by H2 is the island of cooling 
below T = fO'*. Within this island, Aheat is dominated by H-photoionization heating, and 
Acooi is dominated by H2-cooling; the net Atot is negative. Note that the island diminishes for 
a thin shielding slab, while its boundaries expand as the column density of the slab grows. 

The reasons for the inverse relation between the area of the H2-cooling island and 
the total transmitted fiux F = J J^dp are twofold. The equilibrium H2-fraction /H2,eq 
is set by a balance between the two-step coUisional formation process and the direct 
photodissociation rate. Hence, /H2,eq increases with decreasing fiux, and so does the H2- 
cooling. In addition, photoionization heating decreases with decreasing fiux. Roughly, the 
scahngs are Acooi ~ 1/VF and Aheat oc F, so that Acooi/Aheat ~ F~^/^. 

The exponential damping factor in equation 3 reduces the total fiux, and causes the 
boundaries of the cooling-island in Figure I to expand (note that a similar effect occurs 
if the overall normalization, J21, is lowered). However, it should be borne in mind that 
as the ionizing fiux is reduced, the formation rate of H2 is lowered (oc a/F), and chemical 
equilibrium is less likely to be achieved within a Hubble time. In order to address this issue, 
we integrated equation (I) as a function of time and calculated the H2 fraction after a Hubble 
time at 2; = 5, /H2,Hub- Initially, the gas was assumed to be neutral and with no H2, i.e. 
nii-=nii+=niie+ = niie++=nYi^ = 11^^+ = n^- = and nHe/?^H = 0.078. 

Figure 2 shows /H2,Hub and tcooi/^dyn as functions of the shielding slab width in 
units of the Jeans length, for a fixed cloud density and temperature of n = I cm~"^ and 
T = 10^ K. The upper panel demonstrates that as the optical depth of the shielding slab 
increases, the H2-fraction initially grows while maintaining chemical equilibrium. The value 
of /H2,Hub peaks at a slab thickness corresponding to tchem = ^Hub, and declines at yet 
thicker slabs where it is out of equilibrium. The lower panel of Figure 2 illustrates that the 
effectiveness of H2-cooling is maintained at much greater optical depths than allowed by the 
requirement of chemical equilibrium. Indeed, the cooling time-scale remains much shorter 
than the dynamical time even as the slab thickness exceeds the Jeans length by two orders 
of magnitude. Fven if H2 is out of equilibrium, the increase in the optical depth can increase 
the ratio Acooi/Aheat sufficiently to turn heating into cooling that is capable of affecting the 
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gas dynamics. For comparison, Figure 2 shows how heating turns into coohng for three 
different normahzations, J21 = 0.05, O.fO, and 0.f5. 

In conclusion. Figures f and 2 demonstrate that as the radiation field penetrates 
deeper into a gas cloud, it is more likely to cool the gas rather of heat it. There are, 
however, several limitations on the physical applicability of these results. In the rest of this 
section we discuss each of these limitations and demonstrate that they play only a minor 
role in most of the parameter space of the Il2-cooling island. 



3.1. Cooling time vs. dynamical time 



The contour-diagram shown in Figure I provides only the sign of Aheat — Acooi- 
However, the fact that this quantity is negative within the Il2-cooling island does not in 
itself guarantee dynamical consequences. In order for cooling to effect the dynamics of a gas 
cloud, the cooling time needs to be shorter than the dynamical time. The line tcooi = ^dyn is 
shown for the five values of column density, A^h = 0, 10^®'^^'^°'^^ on the density-temperature 
plane in Figure I, where we adopt the following definitions for the cooling time. 



(3/2)nfcBr _ (3/2)/>fcBr//imp 
~|A_AI~IA_AI ' ' 

M-^heat ^^cool M-^heat ^^cool 



and the dynamical time 



It is apparent from Figure 2 that most of the Il2-cooling island, except for its high-density 
high-temperature corner, lies in the region where tcooi < ^dyn- This implies that Il2-cooling 
in most cases is indeed efficient enough to affect the dynamics of the gas. 



3.2. Chemical equilibrium vs. Hubble time 



Figure I assumes chemical equilibrium. If the time-scale to reach chemical 
equilibrium, tchem, is too long, the Il2-fraction could be significantly lower than its equilibrium 
value /H2,eq- Although the chemical reaction-network in equation (I) does not possess a 
single time-constant that characterizes relaxation to chemical equilibrium, the time-scale for 
II2 to reach equilibrium is roughly proportional to the square of the number density of free 
electrons, i.e. oc Ij F. 

We compare tchem and the Hubble time at 2; = 5, tn = 2.0 X 10^^ s (assuming h = 0.5 
and = I). Although the chemical time-scale could also be compared to the dynamical 
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time, this would be less relevant for virialized objects that do not undergo bulk motions in a 
dynamical time. Figure 1 shows (dashed) lines of tchem = ^Hub for three values of the column 
density, A^h = 10^^'^°'^^. For A^h = and 10^®, H2 is in equilibrium every where in Figure 1. In 
the regions where tchem < ^Hub, chemical equilibrium can. in principle, be established. In the 
regions where tchem > ^H, the Il2-fraction reached within a Hubble time /H2,Hub depends on 
the initial abundances and is always lower than the equilibrium abundance /H2,eq if initially 
/h2=0. Nevertheless, we find that the contour-diagram is practically unchanged if /H2,Hub 
is used in place of /H2,eq to calculate Atot- The out-of-equilibrium II2 fraction becomes 
dynamically important, with Atot < and tcooi/^dyn < 1 within a Hubble time, even for the 
lowest column density shown in Figure 1, A^h = 10^^ cm~^. 

3.3. Chemical equilibrium vs. thermal equilibrium 

In realistic circumstances the approach to thermal and chemical equilibrium may be 
coupled. This coupling could affect the interpretation of Figures I and 2. Here, we will 
consider the regime where both the thermal and the chemical time-scales are shorter than 
the dynamical time. We therefore keep the density fixed, and consider the case n = lOOcm""^. 

If the initial H2 abundance is small, then as soon as the background radiation turns 
on the gas will heat. The subsequent thermal evolution of the gas, however, depends on 
the competition between the photo-ionization heating rate and the increase in the H2Cooling 
rate. Both the initial heating rate and the initial H2-formation rate increase with increasing 
incident fiux. However, the heating rate is proportional to the total fiux F, while the 
H2-formation rate grows only as a/F. The H2-formation rate is also sensitive to the gas 
temperature; it increases with temperature up to ~ 3, OOOK but then decreases abruptly at 
yet higher temperatures because of coUisional dissociation. Thus, the relevant question is 
whether photoionization heating raises the gas temperature above ~ 3, OOOK before H2 is 
formed in sufficient quantities to reverse the sign of Aheat — Acooi- As we show below, the 
answer to this question is positive for high incident fiuxes and negative for low fiuxes. 

We solved the coupled time-evolution of the temperature and the relative abundances 
during the approach to thermal and chemical equilibrium. This was done by integrating 
equation (I), starting from T = 10^ K and allowing the temperature and the rate coefficients 
ki^ to evolve at each timestep according to Aheat — Acooi- The initial abundances were chosen 
as in Figure 2. Figure 3 illustrates the sensitivity of the resulting temperature and H2- 
fraction to the incident fiux, characterized by the slab column density A^h- For a higher fiux 
(A^H = 0, upper panel), Aheat remains roughly constant since the ionization fraction is low 
but Acooi is increasing as more H2 is produced. When the temperature reaches ~ 3, 000 K, 
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molecular hydrogen is destroyed. The final equihbrium H2-fraction is correspondingly low, 
/H2,eq = 4.0 X 10~^, as the final equilibrium temperature reaches the high value of 8,440 K. 
The lower panel shows a case when the incident fiux is lower, due to a thick slab A^h = 10^^. 
Here the H2-fraction increases sufficiently rapidly so as not to allow the temperature to reach 
~ 3,000 K. In this case, the final equilibrium H2-fraction is high, /H2,eq = 2.5 X fO~^, and 
the equilibrium temperature is low, 290 K. In fact, the equilibrium temperature in this case 
is below the initial temperature of 1,000 K, implying that a net cooling feedback due to 
molecular hydrogen is possible. 

3.4. Additional Processes 

In addition to the time-scale issues addressed in § 3.1-3.3, we examined three processes 
that could invalidate our conclusions. We found all three processes to be insignificant, and 
therefore only comment on them briefiy below. 

3.5.1 II2 destruction by the Solomon process 

Il2-molecules can be dissociated by photons with energies below I3.6eV by a two- 
step process, as first suggested by Solomon in 1965 (cf. Field et al. 1966) and later shown 
quantitatively by Stecher & Williams (1967). The incident number fiux of photons is likely to 
be higher below the H-ionization threshold of I3.6eV than above it, especially if the radiation 
goes through a slab of neutral hydrogen which absorbs photons only above I3.6eV (cf. eq. [3]). 
Nevertheless, we find that because of its small cross-section cTsw ~ 3 X (I0~^^ — 10~^^) cm^, the 
Solomon process does not effect the contour maps shown in Figure I. As long as there is a 
high-energy component to the incident spectrum, the Solomon process is unable to compete 
with the direct photodissociation of II2 by the high-energy photons. 

3.5.2 Fast electrons 

Il2-molecules could also be dissociated by fast electrons that arise as the product of 
ionization by energetic photons [hp ~ IkeV). In fact, depending on the ionization fraction 
(which determines the drag on the fast electrons) one fast electron can dissociate several II2- 
molecules. Xu and McCray (I99I) give the fractional energy deposition into thermalization, 
excitations, and dissociations, by fast electrons for various energies and ionization fractions. 
Based on these fractions, we conclude that even though as much as 30% of the electron 
energy can go into dissociations, the number of Il2-destructions by fast electrons is still 
small (< 15%) in comparison with the number of direct Il2-photoionizations. Consequently, 
the boundaries in the contour map (Fig. I) would not change even if multiple-dissociations 
by fast electrons were included. 
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3.5.3 Self-absorption of infrared photons 

The H2-cooling could, in principle, be suppressed if the gas cloud were optically- 
thick to its emitted infrared-line photons. H2-cooling occurs mainly due to radiative de- 
excitation from excited rotational and vibrational levels of the electronic ground state. We 
obtained cross-sections for these transitions from the transition probabilities given in Turner 
et al. 1977, and found them to be smaller than ~ 10~^^ cm^. Given these small cross-sections 
and the typically small H2-fraction ( ^ 0.01), we conclude that an exceedingly high H-column 
density, ~ 10^® cm~^, is necessary to absorb the IR cooling photons. 

4. Cosmological Implications 

In summary, we find that a UV irradiation of a cold gas cloud with a density ^ 1 cm~"^ 
can trigger excess H2-cooling on a time-scale shorter than its dynamical time. The effect 
is more pronounced in the inner cores of massive gas clouds which are shielded from 
photoionization heating by the surrounding material (cf. Fig. 1). Contrary to naive 
expectations, the existence of a strong ionizing background at redshifts 2; ~ 2 — 4, could 
have accelerated the collapse and fragmentation of dense objects, rather than increased their 
thermal pressure support. This H2-cooling could have lowered the Jeans-mass threshold 
for fragmentation and led to enhanced star-formation rates in these objects. Since the UV 
background is found to decline to low values at low-redshifts (Kulkarni & Fall 1993), this 
enhancement in the fragmentation efficiency may have been a unique feature of the universe 
at redshifts 2; ~ 2 — 4. 

The H2-cooling effect may also have important implications on the collapse 
of objects, and the duration of the reionization epoch at higher redshifts. The 
cooling island in Figure 1 extends down to a density corresponding to a redshift 
(l + z) 30 (Ob/iV0.025)"^^^ (1 + A/200)"^/^ where A is the overdensity in a collapsing 
region and Ob is the baryonic density parameter. In particular, if the very first objects 
to "light up" triggered the cooling effect discussed here, they could have accelerated 
the formation of additional ionizing sources, and reionization would have been completed 
abruptly. In contrast, if the first sources did not exert this positive feedback (e.g. because 
their UV spectrum was stellar, and did not extend towards the X-ray band to trigger the 
formation of II2), then the photoionization of the IGM must have taken longer, with the 
Jeans mass rising continuously up to 10^~^Mq. These two scenarios differ in the amount of 
small-scale structure formed by the end of the reionization epoch. 
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Table 1: The equilibrium H2-fraction /H2,eq = 2nH2/(?^H + n-u- + n-H+ + 2nH2 + 2nH2 + ) as a 
function of total number density n, temperature T, and column density A^h- For each column 
density, the five rows correspond to five different temperatures 2.0 < log(r/K) < 4.0, the 
eight columns to the number densities — 1.0 < log(n/cm~"^) < 6.0. 
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Fig. 1. — Zero energy exchange contours (Aheat — Acooi = 0) for five different column- 
densities A^H = 0, 10^®'^^'^°'^^ cm~^. The five solid lines in the upper right corner show 
corresponding lines of tcooi = ^dyn (the right-most line corresponds to A^h = 0, the left-most 
line to A^H = 10^^) . Going to the left, the area enclosed by the lines defines increasingly 
larger regions of ineffective cooling (tcooi > ^dyn)- The three dashed curves show lines of 
^chem = ^Hub fof A^H = fO^^ (lowcr curvc) and A^H = 10^°'^^ (upper two curves). Above these 
lines, chemical equilibrium is established in a Hubble time. For A^h = and 10^®, the lines 
are off the plot. 
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Fig. 2. — Effects of sliieiding by a siab of tliickness £ on tlie H2-fraction readied in a Hubble 
time (upper panel), and on the ratio between the cooling time and the dynamical time (lower 
panel). In the lower panel, cooling is denoted by a solid line and heating by a dotted line. 
The gas is assumed to have a density 1 cm~"^ and a temperature T = 10^ K. The horizontal 
axis describes the thickness of the shielding slab in units of the Jeans length (bottom label) 
and in the corresponding column density log A^h (top label). 




Fig. 3. — The approach to equihbrium when thermal and chemical evolution are coupled. 
The vertical axis shows log of the H2-fraction on the left, and log of the temperature in K 
on the right. The heating and cooling functions are arbitrarily normalized. 



